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Abstract. Discrete multiplicative turbulent cascades are described using a formalism involving infinitely 
divisible random measures. This permits to consider the continuous limit of a cascade developed on a 
continuum of scales, and to provide the stochastic equations defining such processes, involving infinitely 
divisible stochastic integrals. Causal evolution laws are also given. This gives the first general stochastic 
equations which generate continuous multifractal measures or processes. 

PACS. 02.50.Ey Stochastic processes - 05.40.Fb Random walks and Levy flights 



1 Introduction 



Multiplicative cascades were first introduced in turbulence 
to model the energy flux in the inertial range. The cas- 
cade formalism was originally introduced as a discrete (in 
scale) procedure, with a fixed (often 2) scale ratio be- 
tween the scale of a structure and that of the daughter 
structures (see [|l|,|[||,^). Discrete cascade models lead 
to discrete scale invariance, characterized by log-periodic 
modulations [^,|6[0,^. On the other hand, a continuous 
symmetry leading to scaling for any scale ratio has been 
proposed, and corresponds to a scale densification of cas- 
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,l^;|l^. Scale densification im- 
plies the use of infinitely divisible (ID) random variables, 
defining cascade models that can be called log-ID ||l0|,^, 
, and have been compared to experimental data in var- 
ious studies ||l|,|l|,[^,|l|,|l|,|3. This then leads to an 
interrogation: discrete cascade models are built using a 
simple recursive multiplicative procedure (see below), but 
what is the continuous limit of this procedure? What is 
the stochastic equation generating continuous multifrac- 
tals? Up to now, the process generated by such contin- 
uous multiplicative cascades has not been explicitely de- 
scribed in the general case; only its statistical moments 
are given |ll3| , or some general relation verified by the pdf 
at different scales This paper aims at detailing how 
ID random measures can be introduced for discrete cas- 
cades; their continuous limit is then given in the form of 
ID stochastic integrals. The stochastic evolution laws that 
generate causal continuous multifractal processes will also 
be provided. This is of direct importance for providing 
estimators of the future state of the process. 



2 Discrete multiplicative cascades 



Soon after Kolmogorov and Obukhov published their log- 
normal proposal for the statistics of the small-scale dis- 
sipation field 1^,^, experimental studies showed that 
the dissipation field had long-range power-law correlations 
[^|2^. This lead Yaglom to propose a random cascade 
model with long-range correlations and small-scale lognor- 
mal statistics . Yaglom's multiplicative cascade model is 
at the basis of most cascade models introduced later to ac- 
count for turbulent intermittency. It is a discrete (in scale) 
model, but most of its properties are shared by continuous 
models. A lognormal pdf is assumed, but this is an unnec- 
essary hypothesis, as is now well recognised. This model 
is multiplicative, nested in a recursive manner. The multi- 
plicative hypothesis generates large fluctuations, and the 
stacking generates long-range correlations, giving spatially 
to these large fluctuations their intermittent character. 

As is classically done (see e.g. Frisch |^), we define 
the cascade yielding a dissipation field e(x) at the smallest 
scale £o, as the product 



(1) 
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of n independent realisations Wi^x of a common, positive 
law. The cascade is developed from the largest scale L 
down to £o — L/A where yl = A" is the total scale ra- 
tio and Ai > 1 is the constant scale ratio between two 
consecutive scales. Generally one assume for convenience 
Ai = 2, but we will later on consider the Ai — > 1 limit cor- 
responding to a continuous cascade |^ . Since all random 
variables are independent, one has the moments of order 
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g > of e: 

n 

1=1 

where K{q) — \ogx-^ < >. For discrete cascades, the 
analytical expression taken by K(q) is only loosely con- 
strained a priori: by conservation K{0) = 0, K{1) = 0, and 
since K{q) is - up to a log Ai factor - the second Laplace 
characteristic function of the random variable log W, it is 
a convex function (see |p8| ). 

To densify the cascade described above, we keep the 
total scale ratio A large but fixed; the continuous limit can 
be obtained by increasing the total step number 1% hence 
Ai = yli/" ^ 1+ (see [|3,|l|,|l|,|9)). Equation @j then 
shows that, in this limit, log e is an ID random variable (see 
psi for ID random variables): continuous cascade models 
are log-ID @|l|,|l|,|lj. This restricts the eligible cascade 
models, since ID laws define a specific family of probability 
distributions. 

Let us also mention one of the main properties of multi- 
plicative cascades: long-range power-law correlations. Fol- 
lowing the development given by Yaglom 0] , one can con- 
sider two points separated by a distance r as having com- 
mon ancestors from steps 1 to p, and separated (hence the 
corresponding random variables are independent) paths 
for steps p to n, where A^ ~ t. Direct calculations then 
provides the classical result for two-points correlations of 
multifractal fieds |pO|,^: 

< €{xYe{x + rf >« yl^(P+'?)r^(p)+^('?)-^<'(P+9) (3) 

for p > and q > 0. Since K{q) is non-linear for multi- 
fractal distributions, the exponent K{p)+K{q)—K{p + q) 
quantifies the long-range power law correlations of multi- 
fractal measures. For p — q = 1, this yields the n — K{2) 
exponent originally given for usual correlations by Yaglom 

The scaling law for the moments (p and the power-law 
correlations (0) are the two signatures of multifractality, 
that are to be recovered by multifractal stochastic models, 
as we define them below. 



3 ID random measures and stochastic 
integrals 

The densification of a multiplicative cascade implies that 
log e is an ID random variable; we now express this den- 
sification in the form of ID stochastic integrals. Since we 
need below to consider moments of order g > of the ID 
random variable F — loge, we consider ID laws for which 
the second Laplace characteristic function ^x{q) — log < 
e"^-^ > converges for a given domain 0. For an ID random 
variable X, one has the general result that Vn integer, 
^^^-^ is still a second characteristic function. This shows 

n 

that a family of ID laws can be defined: two ID laws belong 
to the same family if their second characteristic function 
is proportional. Then each ID family can be character- 
ized by a reference function 3'o('?)- We choose a reference 



function such that !fb(l) = Ij a-nd for an ID random vari- 
able X, we define its scale S{X) as the proportionality 
factor, giving the general identity ^x[q) — S{X) ^o{<l)- 
The scale is a positive real number; it is an additive func- 
tion since for two independent ID random variables X 
and Y of the same family, it is easily checked that we 
\i{weS{X + Y) = S{X) + S{Y). An ID random variable is 
then uniquely characterized by its reference second char- 
acteristic function and its scale (relative to this reference 
function). As examples, ^'0(9) = <f for a Gaussian law, 
and <fb(9) — log2(l -f q) for a Gamma law. 

We then define ID random measures as set functions 
M{A), such that VA, M{A) is an ID random variable, 
with a scale given by S{M{A)) — m(A), m{A) being 
the control measure of M{A). We can easily check that 
M{A) possesses the basic additive property of random 
measures: for two sets A and B with AC] B ~ 0, let us 
note C = A U B. By definition, M(C) is a ID law, and 
its scale verifies: S{M{C)) = m{C) = m{A) + m{B) = 
S{M{A)) + S{M{B)), hence that M(C) = m\a)+M{B). 
This corresponds to the following second characteristic 
function for M{A): 

'FM(A){q)--m{A)^Q{q) (4) 

This expression takes a more familiar form in the ID case 
when A — [0, t] is an interval, and taking M{A) = Y{t) — 
Y{0) = Y{t) where F is a process with independent and 
stationary increments. Then one has the classical result 
^Y{t){(l) = t ^z{<l) where Z is the stationary process given 
by Z{t) = Y{t) - Y{t - 1). In the following we keep the 
random measure notation, which is more general. 

Having introduced an ID random measure M, a stochas- 
tic integral can be built (see e.g. [|3|), as a Stiltjes integral: 

„6 n-l , _ 
/ f{t)M{dt) - hm V/(a + * -) 

°- i=0 

6 — a , ,6 — a,,,, 

M{{a^i ,a+ i-Hl 5 

n n 

As we will see below, the densification of the cascade leads 
to a stochastic integral with /(i) = 1 \lt. In this case, the 
second characteristic function of the integral has a simple 
expression. Let us note / = J^M{dx). By additive prop- 
erty, I is still an ID law of the same family as M, and its 
scale is given by S{I) = S {J.M{dx)) = J^dx ^ rn{A), 
such that we have still Eq. (|) with M (A) = M (dx) 
and m{A) — Jj^dx. Let us note that when / is not identi- 
cally 1, the result is not so simple, since an addition of ID 
random variables belongs to the same family, but not a lin- 
ear combination. One has is this case (if the integral con- 
verges, see Q): ^i{q) = J^^^q {qf{x)) dx, showing that 
in general, 'PJIq) is not proportional to !fb(9)- 

4 Densification of the cascade and stochastic 
equations 

Let us introduce A a variable scale ratio, verifying 1 < 
A < yl, where A is the fixed largest scale ratio. We also 
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introduce R = log A and r — log A. The elementary scale 
ratio introduced above writes now Ai — A^^" = e^/". The 
discrete cascade corresponds to introducing a stochastic 



kernel M and intervals Ap and Bp such that 



r{x)^logeix) = {Ap,Bp{x)) 

p=0 



(6) 



where here M(A, B) is a random variable depending only 
on m{A), giving ^M(A,B){q) = 'm{A)^a{q)- The intervals 
Ap and Bp, responsible for the cascading parent/children 
structure, are built in the following way: the width of Ap is 

linear in r, giving Ap = '^^^^^^' ]. The intervals Bp{x) 
are centered in x and of width proportional to A^ = gpR/n-^ 
giving Bp{x) = [x - fe^^/", x + §6^^/"] where t = L/A 
is the resolution. The densification, corresponding to n — > 
oo, transforms then Eq. into a stochastic integral and 
using Eq. (|^) (with /=1), we obtain finally that: 



e/i(a;) — A exp 



Ml ^,i?,/o(x) 



(7) 



where c > is a parameter, Iq [x) is the interval of length r 
centered in x, and D\ is the dilatation operator of factor A. 
At a given position a;, the stochastic integral corresponds 
to a kernel visiting a conical structure, as represented in 
Fig. 2. This expression can be generalized to d-dimensional 
domains, and also to anisotropic scaling symmetries. 

It can be easily verified that this stochastic equation 
generates a multifractal field. Indeed, the moments are 
scaling as < e\ >= A-^^'^^ with 



Kiq) = ci%iq)-q) 



(8) 



Moreover the two-points statistics can also be recovered: 
as was done above for Eq. (^, the correlation < (eA(x))^ 
(e/i(a; + y))'^ > involves two integrals which have no inter- 
section (and thus are independent random variables) for 
A < Ao , where Xq = = ^. Aq is the scale ratio of transi- 
tion, as shown in Fig. 3, and for Aq < A < yl, the random 
variables corresponding to the two stochastic integrals are 
no more independent. After some calculations, this leads 
to the same expression as Eq. (||) for discrete cascades. 

We also give the expression for causal cascades, where 
the position is time and the past does not depend on the 
future. This case is of particular importance for predic- 
tion of multifractal times series. In this case, one can 
readily modify the intervals Bp{t) by taking an interval 
of the same length as before but preceding t: Bp{t) = [t — 
T exp{pR/n),t]. This gives the following causal stochastic 
evolution law for continuous multifractals: 



eA{t) = A-'e^p 1^ M ^ ^ 



f cdX 



[t - rA, t] 



(9) 



Let us finally consider an irnportant family, correspond- 
ing to logstable multifractals |||,^, including the lognor- 
mal case. Stable laws are ID, and possess a stronger prop- 
erty corresponding to stability, which can be written here 



M{kA) = fci/"M(A) for fc > 0, where < a < 2 is 
the Levy index; a — 2 for the Gaussian case 1 33 , 3^ . We 
have here ^o{q) — q"; we note that when a < 2, the sec- 
ond Laplace characteristic function is defined for positive 
moments only for asymmetric laws for which hyperbolic 
pdf corresponds to negative fluctuations (Pr{—X > x) ~ 
x~"), whereas positive fluctuations have an exponential 
decay Then, by splitting Eq. (||) into two inte- 

grals, corresponding to backward and forward domains, 
and introducing the change of variables u = x — and 
V = x-\-^X respectively, one obtains (introducing the Levy 
measure La{du) — M{du, [u,x\)) a stable stochastic inte- 
gral: 



eA{x) = A- 



' exp / \u- 

I A{x) 



(10) 



where A{x) = [x - X/2, x - t/2]U[x + t/2, x + X/2]. This 
equation corresponds to the exponential of a fractional 
integration (over a limited domain) of order 1 — 1/a of a 
Levy-stable noise. This expression was already given 
by phenomenological arguments, and is here derived 
as a direct consequence of the densification. When the 
position is time, we obtain the following causal stochastic 
evolution equation for a logstable multifractal generated 
with a fixed scale ratio A = T/t: 

ft-r 

CAit) = A-'^ exp {t-u)-^'"dLc,{cu) (11) 

Jt-T 

This can be directly used in numerical simulations. For 
lognormal multifractals La is replaced by the Wiener mea- 
sure W. 



5 Conclusion 

We now discuss the new results provided by our approach, 
compared to related papers. She and Waymire |13| have 
given a general expression for K{q) for continuous cas- 
cades, using the canonical Levy-Khinchine representation 
for the second characteristic function of ID laws. This 
could of course also be provided here. On the other hand, 
we have pushed further the analysis, since the process it- 
self was not studied in Castaing and collaborators 
have studied the convolution properties of the probability 
density of velocity increments under a change of scale, for 
continuous cascades [^,|l^,^. This property is recovered 
here (it can be obtained from Eq. (0)), and we also pro- 
vide the stochastic equation for the process itself. Other 
studies have provided Fokker-Planck or Langevin [ p6| 
equations for the cascade process in the scale-ratio space. 
These equations apply only to lognormal cascades. These 
studies consider a fixed spatial (or temporal) position and 
provide stochastic equations for the cascade process devel- 
oping at this particular position. This framework is then 
different from ours by the fact that (i) the position is not 
taken into account, whereas we included the position in 
our equations, and (ii) we provided equations for the gen- 
eral log-ID case, and not only for the lognormal case. 



4 



Frangois Schmitt, David Marsan: Stochastic equations generating continuous multiplicative cascades 



We have obtained continuous niultifractals as the ex- 
ponential of a stochastic integral. This formalisation is of 
great interest for theoretical studies since it provides the 
first general stochastic equations for continuous multifrac- 
tals. These equations can be generalized to anisotropic sit- 
uations, and their dynamical properties (i.e. error growth 
and the corresponding predictability limit, or return times) 
are now open to theoretical studies, whereas before such 
studies were possible only through numerical simulations. 
We have also provided the general evolution equation for 
causal continuous multifractal processes. We have shown 
how these equations simplify for log-stable and lognormal 
multifractals. 

These equations can be used for numerical simulations 
of continuous multifractals. In the general log-ID case, Eq. 

can be used as follows: a 2D ID noise must be gener- 
ated, in the r {r — log A is the logarithm of the scale 
ratio) - position plane. Then for each position, a numeri- 
cal path-integration is performed in this plane, as given by 
this equation. In the stable case the procedure is simpler, 
since it is enough to simulate a ID stable (or Gaussian) 
noise and to proceed to a fractional integration over this 
noise, as given by Eq. (p^. 
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